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ABSTRACT 

Within the context of constraining an expansion of the dark energy equation of state 
w{z) we show that the eigendecomposition of Fisher matrices is sensitive to both the 
maximum order of the expansion and the basis set choice. We investigate the Fisher 
matrix formaUsm in the case that a particular function is expanded in some basis set. 
As an example we show results for an all sky weak lensing tomographic experiment. 
We show that the set of eigenfunctions is not unique and that the best constrained 
functions are only reproduced accurately at very higher order > 100, a tophat basis 
set requires an even higher order. We show that the common approach used for finding 
the marginalised eigenfunction errors is sensitive to the choice of non-u'(z) parameters 
and priors. The eigendecomposition of Fisher matrices is a potentially useful tool that 
can be used to determine the predicted accuracy with which an experiment could 
constrain 'w{z). It also allows for the reconstruction of the redshift sensitivity of the 
experiment to changes in w(z). However the technique is sensitive to both the order 
and the basis set choice. Publicly available code is available as part of iCosmo at 
http : //www. icosmo . org. 
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1 INTRODUCTION 



] The Fisher matrix methodology (Tegmark, Taylor & Heav- 
e ens, 1997; Jungman et al., 1996; Fisher, 1935) is a statistical 

■ tool that has been used with some success in predicting the 
, ability of future experiments to constrain particular param- 

■ eters of interest. In cosmology Fisher matrices have gained 
' some importance in predicting the potential outcome of ex- 
, periments, in particular dark energy surveys, on which a 
' large amount of resources may be spent. It is therefore of 

paramount importance that the way in which Fisher matri- 
ces are used should be understood and that any results that 
depend on this methodology should be robust and reliable. 

In this article we will outline how the techniques of de- 
composing a general matrix have a special interpretation 
when used within the Fisher matrix framework. This will be 
done within the context of attempting to predict constraints 
for an extended parameter set. We will use the specific ex- 
ample of attempting to constrain dark energy equation of 
state w{z) parameters (Albrecht & Bernstein, 2007; Tang et 



al., 2008; Simpson & Bridle, 2006; Crittenden & Pogosian, 
2005; de Putter & binder, 2008; Dick et al., 2006; Huterer & 
Cooray, 2005; Huterer & Starkman, 2003; Hu, 2002; Knox 
et al., 2005; Ishak et al. 2006; Zhan et al., 2009a; Zhan et 
al., 2009b; Cunha et al., 2009; Joudaki et al., 2009; Sarkar 
et al, 2008), most recently Albrecht et al. (2009) defined a 
new dark energy Figure of Merit using a binned 36 param- 
eter model. 

We will show that the eigenfunctions obtained by diagonal- 
ising a Fisher matrix are dependent of the basis set used in 
the construction of the Fisher matrix, and that the eigen- 
functions only tend to agreement in the limit of a very large 
order of expansion. We will also show that the usual method 
of marginalising over extra parameters - constructing the 
Schur complement (Zhang, 2005) of the total Fisher matrix 
- is dependent on the choice of non-w^z) parameterisation 
and priors when also performing an eigendecomposition. 

In Section 2 we will introduce the methodology, in Section 
3 we will apply the approach to weak lensing tomographic 
survey prediction, we present conclusions in Section 4. 
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2 METHODOLOGY 

For a set of parameters 9 the Fisher matrix allows for the 
prediction of parameter errors given a specific experimental 
design and method for extracting parameters. In the case of 
Gaussian-distributed data where we assume that the error 
on the signal is not a function of parameter values gc 
ac{0) we can take the covariance of the estimated values of 
the parameters; 

cov[0„0,] = {(0, - - {0.»> = (1) 

where the Fisher matrix is defined by (Tegmark, Taylor & 
Heavens, 1997; Jungman et al., 1996; Fisher, 1935) 

where 5* is the signal. The predicted marginal errors on the 
parameters are given by ASi = ^ {F~^)ii, this is the mini- 
mum marginal error that one can expect for the method and 
experimental design considered (due to the lower bound in 
the Cramer-Rao inequality). 

If we extend the parameter set to </> = (6, ip) then the Fisher 
matrix is extended to include these extra parameters 

and the predicted marginal errors on the parameters A9i = 

ij (F~i)^'^ are increased due to degeneracies between the 
original parameters and the new parameters. 

For dark energy cosmology prediction the parameter set is 
commonly divided into those that contain information on 
the dark energy equation of state w{z) and those that do 
not. The usual approach taken (Albrecht et al., 2009; Al- 
brecht & Bernstein, 2007; Tang et al., 2008; Simpson & Bri- 
dle, 2006; Crittenden & Pogosian, 2005; de Putter & Linder, 
2008; Dick et al., 2006; Huterer & Cooray, 2005; Huterer & 
Starkman, 2003; Hu, 2002; Knox et al., 2005; Ishak et al. 
2006; Zhan et al., 2009a; Zhan et al., 2009b; Cunha et al., 
2009; Joudaki et al., 2009; Sarkar et al., 2008) is to "bin" 
w{z) in redshift and to assign an amplitude to the value 
of w{z) in each bin. However w{z) is not a directly observ- 
able quantity but can only be inferred through other integral 
relations, most simply the comoving distance or Hubble pa- 
rameter 

Diz) ^ 



Ho 



Az 



(n„(l + Z)3 + n^ee^SS d.'[l+»{z')l/(l + .'))i 



(4) 



H{z) ^ Ho{fl,r,{l + zf + n^.e' 



3/o"dz'[l+»(z')]/{l+z')Nl 



(5) 



So any binned expansion of 'w{z) to include extra parameters 
cannot be justified by arguing that the data itself is binned 
in redshift. 

The choice of how to expand the function 'w{z) is thus a 
purely theoretical one. In general we will consider ex- 
panded in some complete basis set 



where an are the basis coefficients and (fin are the basis func- 
tions. Since we do not know what dark energy equation state 
is apt for the physical Universe any chosen basis set is as 
valid as any other. Binning is a special case where the basis 
functions are tophat functions in redshift 

0^"'"-(^)=H(z-^)-H(z + ^), (7) 

where Az is the bin width and TC is the Heaviside (step) 
function. 

Note that interpolations between delta function values at 
discrete redshifts, linearly or with some spline function, are 
different basis sets to binning. We clarify three broad mean- 
ings of binning that occur in the literature 

• Tophat basis set : functions are not continuous, all 
derivatives are ill-defined. 

• Linear interpolation : functions are continuous, I''' 
derivatives are discontinuous, 2"^^ and higher order deriva- 
tives are ill-defined. 

• Spline (e.g. quadratic spline) : functions are continuous, 
1°' derivatives are continuous, 2""* derivatives are discontin- 
uous, 3'^'* derivatives and higher are ill-defined. 

Of these three basis functions, only tophat form an orthog- 
onal basis set; though all form a complete basis set. We 
include interpolation in Section 3 for completeness, and be- 
cause it is commonly used in the literature. 

The most general condition that we could apply to the ex- 
pected w{z) is that it is continuous (features could be very 
sharp in redshift but not physically discontinuous - the 
sound speed of dark energy is generally expected to be < c)^ . 
We note that binning is not discrete differentiable since at 
the bin boundary the gradient is infinite. 

We will now investigate how the Fisher matrix should be 
treated when such an expanded basis set is used. 



2.1 Basis set decomposition 



The general framework with which we are now presented 
consists of a Fisher matrix 



F : 



pee pew{4>) 



(8) 



where are the cosmological parameters not associated with 
wi^z). 'w{(fi) are the set of parameters that describe w{z) 
where we have expanded w{z) using a general complete or- 
thogonal basis set with basis functions (p. The sub-Fisher 
matrix for w{z) is calculated by 

^ a(a;)2 dai daj ' ^ ' 



-,W(4,)W(4,) 



1 -I- W{z) = ^ an4>n{z), 



(6) 



t There are some tachyonic dark energy models (e.g. Bagla et 
al., 2003) but choosing a basis set based on these specific models 
would be optimistic. 
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maximum order in Cheh maximum order in Lege maximum order in TopH 

Figure 1. An example of the residual that can occur when reconstructing a given basis sets functions with another basis set. On the y axes 
are the order of the basis function to be reconstructed, on the x axes are the maximum order of the series that is attempting to reconstruct 
the functions. The gray scale shows the maximum residual (black low residual, whit large residual), each function is investigated over 
the range x G [—1,1]. The left plot shows the order of Chebyshev polynomials needed to reconstruct Legendre polynomials. The middle 
plot shows the order of Legendre polynomials needed to reconstruct Chebyshev polynomials (the inverse operation to the left plot). The 
right hand plot shows the order of the tophat basis set needed to reconstruct Chebyshev polynomials. The red (gray vertical) lines are 
at AT = 36 (Albrect et al., 2009) to guide the reader. 



where (7{x) is the expected marginal error on the signal S 
and the errors on the coefficients of the expansion (equation 

6) yCP^^^^^^)^ = Aoi (not marginalising over non- 
w{z) parameters G). 

A basis set is orthogonal if it satisfies the orthogonality re- 
lationship 



(t)m{z)ct>niz)M{z)dz = Cm(5„„, 



(10) 



over a range R, where M{z) is a weighting function, c,„ are 
constants and is the Kronecker delta. The coefficients 
needed to construct an arbitrary function f{z) using the 
basis set (p are given by 



(11) 



an = — f(z)(l>„{z)M{z)dz. 



By performing this transformation the non-w(z) errors can 
also affected. There a three possibilities; either F is the full 
w and non-ui Fisher matrix in which case the non-ui parame- 
ters are clearly affected, or F is the ro-only Fisher matrix, or 
F can be the marginalised lu-only Fisher matrix. We discuss 
these later two options in Section 2.3. 

The new Fisher matrix A** is a diagonal matrix. The Q ma- 
trices are called the eigenmatrices (or vectors) which are a 
special kind of Jacobian matrix in which the original matrix 
is rotated such that is diagonalised (see Appendix A). In 
performing this operation we have constructed a new set of 
basis functions, e, that are linear combinations of the orig- 
inal basis set, (f). w{z) can now be reconstructed using two 
equivalent forms 



1 + w{z) = ^ an(l>n{z) = ^ 6„e„(z) 



(13) 



In general we can construct the Fisher matrix (equation 
8) using an orthogonal basis set for 'w{z) (orthogonal with 
respect to a weight function M{z)) but the sub-matrix 
pW(tt>)w{ii>) diagonal since w(z) can, at a min- 

imum, only be observed through integral relations^. 

A common approach is to look for the eigenfunctions of 
w{z). The motivation for this is that the eigenfunctions are 
thought to form an orthogonal basis set. These functions are 
found by rotating the Fisher matrix such that it is diago- 
nalised using the transformation 



^0 
F 



Q^FQ 



(12) 



where Q is an orthogonal matrix, Q Q — I. 



where N is the order of the sub-Fisher matrix p'^W^W _ 
we have introduced this maximum order since Fisher matri- 
ces have, by definition, a finite order. The new functions, e, 
can be constructed using the eigenmatrix 



where (Appendix A) 



(14) 



(15) 



The new Fisher matrix has diagonal elements that are re- 
lated to the errors on the new functions' coefficients 

dS dS ,K 1 



A-T 



T — 



a{x)^ dbi dbj 



{Ah 



(16) 



§ Most articles in the literature (e.g. Albrecht et al., 2009; 
Huterer & Starkman, 2003; Crittenden & Pogosian, 2005) neglect 
the weight function; we note that the tophat basis set is peculiar 
in that it is orthogonal with respect to any weight function. 



In performing such a rotation we have effectively created 
a new basis set that is orthogonal with respect to a new 
weight function (equation 10) that takes into account the 
covariance in the Fisher matrix. 
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We note that an operation of the form F — QAQ^ (or the 
inverse FQ — A) leads to a unique Q if A and F are fixed. 
However in general a further operator P can be applied to 
the new diagonal matrix PA = P{Q^FQ) to map a matrix 
F to any diagonal of the same dimension, where P is also 
a diagonal. P acts like a 'stretching' (or compression), not 
a rotating, operator. P does not correspond to a change 
in parameters but a change in the errors on an eigenbasis 
parameter set. We do not consider this in the remainder of 
the article. 

We also note that in general there exists operators that can 
rotate from a large matrix to a smaller one. For example an 
n X n diagonal matrix A''''^^" can be mapped to a smaller 
mxm diagonal matrix A°™^" where m < n via the operator 
h-tT^^ = AkiA^^^^^Aj^. where Aki is an n x m matrix. So if 
the size of the matrices are not specified then the operation 
F = A(QAQ^)A^ = RAR'^ (where R = AQ) can map a 
matrix to a diagonal of smaller dimension. This is equivalent 
to mapping the basis set expansion from a particular order 
to one of a smaller maximum order. Note that this is in 
general a one-way operation since information is lost in the 
minimisation. In the remainder of this article we will only 
consider mapping between matrices of the same dimension. 



2.2 The non-uniqueness of the eigenfunctions 

It is often assumed that an eigenfunction decomposition 
yields a set of functions that are unique (if they are not 
degenerate; Albrecht et al., 2009; Crittenden & Pogosian, 
2005; Huterer & Starkman, 2003). This statement comes 
from matrix theory; the diagonalisation of a particular ma- 
trix can only be done one way. However this does not mean 
that the eigenfunctions of w{z) from a finite Fisher matrix 
are unique, since we can start with different basis sets and 
this choice is arbitrary. This break-down is due to the finite 
number of coefficients - so is really a convergence statement. 

The general minimum assumption case we are dealing with 
is one in which we consider two Fisher matrices calculated 
by expanding w{z) in two different basis sets (j) and i/'- 

We can only transform one Fisher matrix F'*' F'^ to 
another by changing the basis sets if all the basis set (j> 
functions 1 ^ N can be described using only the functions 
1 —> N from the set tp, and vice versa. For example to rotate 
from basis (j) to tj) then 

AT 

iJi{x)=Y,J^J'PAx) V l<i<iV (17) 
J = l 

needs to be true. However for any conceivable complete basis 
set for finite A'^ this sum will always have some residual. This 
residual Ri{x) can be quantified by 

N 

n = JV + l ■' R ' 



The notation we use throughout is that Q are explicitly 
eigenmatrices, J are explicitly general Jacobian transfor- 
mations, 7? are residuals and F are Fisher matrices. Note 
however that an eigenmatrix Q is a special kind of Jacobian 
J where the result of the transformation is a diagonal. 

In Appendix B we show that if the residual is non-zero then 
the eigenvalues are not equal. 

If the residual between any two basis sets is small, in the 
limit of large A*', then the eigenfunctions will agree to some 
accuracy, and one would expect the best constrained eigen- 
functions from any basis set to converge. However for a finite 
A'' the eigenfunctions will only all agree ffthe Af* basis func- 
tion in one set can be reproduced using only the functions 
1 ^ A'' in other basis set. 

Figure 1 shows the residuals (equation 18) between a vari- 
ety of basis sets as a function of order (in equation 18 we 
replaced oo with 10^ for which we find numerical conver- 
gence of the results in this Figure since 50 < 10''). This plot 
shows that in order to have a residual of < 0.01 when recon- 
structing Legendre polynomials with Chebyshev the order 
of the Chebyshev needs to be at least as high as the order 
of the basis function that is being reconstructed. However 
in the opposite case the low order Chebyshev polynomials 
are difficult to reproduce with Legendre, and the Cheby- 
shev polynomials are difficult to reproduce with the tophat 
basis set (binning). In fact Cheyshev functions are more dif- 
ficult to reproduce (especially at low order) because they are 
bounded in the y direction as well as the x direction. 

If one chooses a maximum order of AT = 36 (Albrecht et al., 
2009) for example then there would be > 0.01 residual in 
reconstructing Chebyshev, Legendre and tophat basis sets 
from each other. As a result the reconstructed eigen func- 
tions would not agree. We will investigate this further in 
Section 3. 

To summarise this section 



• The eigenfunctions created using two different basis sets 
and a finite order will only agree if the residual between the 
reconstruction of the basis functions is zero. 

• For a finite order there will always exist some eigen- 
functions from a given basis set that cannot be reproduced 
using another. So for finite order the set of eigenfunctions is 
non- unique. 

• For a large finite order we expect the best constrained 
eigen functions to tend to agreement, and for this agreement 
to improve as the order is increase d1. 



^ The full proof of this involves Weyl's matrix inequality and is 
beyond the scope of this article, we demonstrate convergence in 
Section 3. 
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2.3 Marginalisation 

Another difficulty in constructing robust eigenfunctions is 
deciding what to do with the non-i«(2;) parameters. A com- 
mon approach is to marginalise over these other parameters 
and then perform the eigenfunction decomposition. This in 
effect mixes all of these parameters with the equation of 
state expansion. 

Referring to equation (8) if we want to know the errors on 
expanded w{z) parameters then we can marginalise over the 
other parameters by constructing the Schur complement 
(Zhang, 2005) of the sub-matrix 



W(4>)W{4,) 
IS 



= F 



pW{4,)e i^poeyip. 



(19) 



The new Fisher matrix ^^'^ f*^) then give errors on 
the w{z) parameters already marginalised over the extra pa- 
rameters. 



We can now diagonalise T?^^'*)''''''^' like 



diag(F,'^(«'^(«) = 



^Qinixcd ^ T ^10 {(^>)iy (0) Qinixcd 



- (g 

= A" 



mixcdNT^ei(;((^>) 



where A™"""* is a diagonal matrix and Q™"""* is the eigen- 
matrix of the Schur complement. Note that it is the sum of 
the terms on the right hand side that form a diagonal not 
necessarily the individual terms themselves. 

Alternatively we can diagonahse pW{,p)w{,j>) .^yj^hout 
marginalising so that (^Q^f pW^w^Qn^ ^ ^» ^-^^^^ ^» 
is diagonal. In this case we can construct the new full Fisher 
matrix in equation (3) for which the TT'^^W'^t'^) component 
is now diagonal 

pee f^Qw^ p0w{4,) 



F = 



(21) 



the cross terms need to be modified so that the correlations 
between 9 and w[cj>) axe correct. We can now write the new 
marginalised errors on the new diagonalised w{(j)) vectors 
like equation (19) so that 



2S 



- ^Qw-^T pew(4,) ^pOOyi pW(4,)e Qw ^^2) 

Note that the {Q^)'^ p'^'-'t')^ effectively acts like rotation on 
{P^O)-^ and that since the w[(j)) are independent param- 
eters to the there is no reason to generally expect that 
the second term is diagonal - why should (^q^)T p^We 
the eigenmatrix of (F^^)~^ when the are independent of 
the w{4,)l Hence the ^ diag(f,^(«'^(«) in gen- 

eral and the full eigenfunctions are not equal to the w[(j)) 
eigenfunctions Q™"''^'^ 7^ Q™, comparing equations (20) and 
(22). 

If equation (20) is used to create eigenfunctions then the 
reconstructed functions will contain the eigenfunctions of 
the ■w{z) parameters mixed with the eigenfunctions from 
the non-u;(2) parameters. 

Our concern is that if the eigenfunctions are mixed then any 



10* 
10° 

10° 

10-' 





- / 
// 

J Fixed Errors 


Fixed Errors 




Marginalised Errora 


Marginalised Errora. 



20 40 60 80 20 40 60 80 100 
N N 

Figure 2. The lines show the predicted fixed errors on the basis 
set coefficients for the original basis set (^j^ and for the rotated, 
eigenfunctions, A]v for a maximum order of 100. We have ordered 
the errors in increasing order of magnitude for clarity, so that 
the function in the top panel does not corresponding to the 
basis function, but is the Af"' best constrained. The top two 
panels show the error for iu(z) not marginalised over any other 
cosmological parameters. The bottom two panels show the errors 
(20) including marginalisation over the other ■aoa-w(z) parameters us- 
ing the method outlined in Section 2.3 equation (22). The key to 
the line styles is given in Table 1 - here we show the Fourier, 
Chebyshev and tophat basis sets. 



10° 
lO'' 
10* 
10'' 
10^ 

18= 

10' 
lO* 
10^ 
lO'^ 
10' 



Ordc3r=10 


Ordcr=30 

— ' 


0rder=60 


Order=100 



Figure 3. The 5 largest eigenvalues as a function of the maximum 
order of expansion. Each line represents a different basis set, the 
key to the line styles is given in Table 1. 



statements on the redshift sensitivity to wiz') are dependant 
on the noxi-w{z) parameterisation and the true w(z)-only 
sensitivity is masked. For example if flatness was assumed, 
if a spectral index is included or its running, or if massive 
neutrinos were included (to name a few cases) then conclu- 
sions made using the mixed eigenfunctions will be different. 

What we propose as an alternative is to find the i!;(z)-only 
eigenfunctions - which are robust to the non-iLi(2:) parame- 
terisation - and to include the marginalisation in a consis- 
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N = 2 N = 'l 
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Figure 4. The eigenvalues as a function of the maximum order in 
the basis set expansion. The key to the different basis sets lines 
is given in Table 1. The upper panels show the eigenvalues for 
the second and fourth best constrained eigen functions and not 
marginalising over the non-io(z) parameters. The lower panels 
show the eigenvalues for the second and fourth best constrained 
eigenfunction marginalising over the non-w{z) parameters using 
the method outlined in Section 2.3 equation (22). 



tent way to find the marginalised errors on the ■u;(z)-only 
eigenfunctions 

• Find the eigenfunctions of w{<l)) by diagonalising F'''. 

• Calculate the new full matrix like equation (21) . 

• Find the marginal errors on the iy(0) eigenfunctions by 
marginalising over the other parameters like equation (22). 

We note that using this approach the errors on the eigen- 
function coefficients will be correlated after marginalising 
over the non-'w{z) parameters. But that these marginalised 
eigenfunction errors will be robust to the non-ui(z) choices 
made. 



priors to any parameters. In Albrecht et al. 2009 they sug- 
gest adding a prior of unity to all w(z) binned parameters. 
We note that this does not have an effect on the eigenfunc- 
tions and in the case of survey optimisation and comparison 
adds a common floor to all scenarios and does not affect the 
relative merit of any survey/method. 



3 APPLICATION 

We will now investigate the effect of choosing different basis 
sets for the expansion of w{z). We will use the example of 
weak lensing tomography (described in Amara & Refregier, 
2007) where a photometric survey is split into redshift bins 
and the auto and cross correlation power spectrum of the 
shear fields are used to infer cosmological parameter likeli- 
hoods. We will present predictions for a Euclid/DUNE-like 
survey (Refregier, 2006) - a 20, 000 square degree photomet- 
ric survey in 5 bands with a number density of 35 galaxies 
per square arcminute and a median redshift of 2 = 1, we use 
10 tomographic bins in the range < z < 5.0. We have ex- 
tended the publicly available iCosmo (Refregier et al., 2008; 
Kitching et al., 2009) package to include basis set expan- 
sion of ^1(2), these additions will be available in version 1.2 
(http : // icosmo.pbworks . com). 

We will investigate 6 different basis sets summarised in Table 
1. We expand these basis sets to a maximum order such that 

Order 

l+w{z) = J2 0'N(t>N{z) (23) 
JV=1 

where the fiducial cosmology is ACDM i.e. the fiducial values 
are ajv = for all basis sets. In performing such a calculation 
the fiducial 'w{z) function must be one that all the basis sets 
can reproduce. 



3.1 Parameter errors 



2.4 Priors 

One can add a prior before the eigendecomposition or af- 
terwards. Depending on the nature of the prior this can 
affect the eigenfunctions themselves. Adding a general non- 
diagonal prior before diagonalisation will act to modify the 
eigenfunctions in a similar way as described in Section 2.3. 

If we add a diagonal prior P (no correlation between the 
prior errors) before diagonalisation then in general this 
will affect the eigenfunctions since ff Qfg^ = Af then 
QFQ^ + QPQ^ 7^ diagonal. There is a special case of adding 
a unity matrix as the prior, in this case the eigenfunctions 
are unaffected since QFQ'^ + QIQ'^ = QFQ"^ + 1 = Af + I, 
and adding the prior before or after the diagonalisation has 
the same effect. 

Throughout the remainder of this article we will not add any 



Figure 2 shows the expected 'w{z) conditional errors on the 
original basis function coefficients and on the eigenfunctions. 
It can be seen that by rotating into an eigenbasis set the 
distribution of errors is skewed to produce more well defined 
functions as well as more poorly constrained functions in 
comparison to the non-rotated set. 

It can be seen from the left panels of Figure 2 that the best 
constrained eigenfunctions have a similar error (the variance 
between the lines is small) but for the poorest constrained 
functions the variance between the basis sets is larger. This is 
because the most well constrained eigenfunctions are similar 
for each basis set whereas the poorest constrained have very 
different functional forms. As shown in Section 2.2 the set 
of eigenfunctions is not unique. 

We also show the errors marginalised over the non-«;(2) pa- 
rameters. For the eigenvalue case this is done using equa- 
tion (22) where the eigenfunctions are the u'(2)-only eigen- 
functions. In all cases the best constrained eigenfunctions 
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Basis 


Functions 


Range R 


Weight M(x) Plot key 


Fourier 


cos(n/r) & sin(na::) 


[— 7r,7r] 


1 


Chebyshev 


Tn{x) = cos(nacos(x)) 


[-1,1] 


(1 -x2)-3 


Legendre 




[-1,1] 


1 


Laguerre 




[0,oo] 


e~^ _ . _ . _ 


tophat 


nix-^)-H{x + ^) 


[—00,00] 


1 


interpolation 


mx + c 


[—00,00] 


1* 



Table 1. A summary of the basis sets investigated. The sets are complete over a range R and orthogonal with respect to the weight 
function M{z). Tophat is equivalent to binning. Interpolation means a linear interpolation between delta functions where m = 
and c = yi — mx\ for the nearest two delta functions x\ < x < X2 to x. * interpolation is not an orthogonal basis set, we include it here 
since this is a commonly adopted expansion. 




K = 4 Order=10 : 


N=4 0rder=30 


N=4 0rder=60 


N=4 Ordcr=100 







0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 2.0 



0,0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 2.0 



Figure 5. The second and fourth best constrained cigcnfunctions for the weak Icnsing tomographic survey considered. We show how the 
functions change as the maximum order in the expansion is varied. The thin lines are for each individual basis set - the key is given in 
Table 1. The thick black solid line shows the mean function over all basis sets at each order. 
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Figure 6. The difference between the eigenfunctions for each basis set and the mean eigenfunction (over all basis sets) for each order. 
We show this residual for the second and fourth best constrained eigenfunctions for the weak lensing tomographic survey considered. 
Each panel matches to the corresponding panel in Figure 5 and the key is given in Table 1. 
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are those that are also the most degenerate with the other 
cosmological parameters, and marginahsing over these can 
increase the error on the eigenfunctions by an order of mag- 
nitude. 

Figure 3 shows how the largest 5 eigenvalues change as a 
function of the maximum order in the expansion. It can be 
seen that a stable regime is found in which the values do 
not change substantially as the order is increased. We ex- 
pand this in Figure 4 where we show how the eigenvalues for 
the second and fourth eigenfunctions change as a function 
of maximum order. The Fourier and Chebyshev values agree 
at higher order, the other basis sets - particularly interpo- 
lations and tophat - are very slow to converge. Even for a 
maximum order of A*' = 100 there remains a large variance 
between the basis sets and that this variance increases as 
the eigenvalues decrease. The eigenvalues are related to the 
eigenfunction errors we show these errors in 

Figure 2. The best constrained A'' = eigenfunction has a 
similar error (for tophat, Chebyshev and Fourier basis sets) 
but by A*' > 2 the errors can vary by a factor > 10. 

There is a large range of eigenvalues between the basis sets, 
the most notable outliers being the tophat and interpola- 
tion basis sets. This is related to the fact that the residual, 
even at a maximum order of 100, between the reconstructed 
eigenfunctions and the true eigenfunctions is significant. We 
investigate this further in Section 3.2. When we marginalise 
over the non-w(2) parameters the variance between the ba- 
sis sets remains but each basis set's errors are effected in 
different ways since the basis functions are degenerate with 
the non-w{z) parameters to different degrees. 

We have performed a numerical test that checks that the 
eigendecomposition of the Fisher matrix is working cor- 
rectly. First we find the eigenfunctions by diagonalising the 
w{z) sub-Fisher matrix. We then use the eigenfunctions as 
a new basis set. We create a new Fisher matrix which takes 
derivatives with respect to the coefficients of the new (eigen- 
basis set). This new Fisher matrix should be diagonal and 
should have diagonal elements that are equal to the origi- 
nal Fisher matrice's eigenvalues. We find that the code used 
in this article successfully passes this numerical check. The 
code used is available as part of iCosmo vl.2 and later. 



3.2 Eigenfunctions 

The eigenfunctions are reconstructed from the eigenmatrix 
Q using equation (14). Figure 5 shows the second and fourth 
best constrained eigenfunctions for different maximum or- 
ders in the expansion for each basis set. In the absence of 
any way to reconstruct the 'true' eigenmatrix we also show 
the mean eigenfunction averaged over basis sets for each or- 
der - in the limit of a high order we expect the eigenfunctions 
to converge. 
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Figure 7. The variance between the eigenfunctions, recon- 
structed using a maximum order in the basis set expansion, and 
the mean eigen function using an order of 100. The key to the dif- 
ferent basis sets lines is given in Table 1. The upper panel shows 
the variance for the second best constrained eigenfunctions and 
the lower panel show the variance for the fourth best constrained 
eigen function. 



nificantly from the mean. We show the difference between 
each basis sets eigenfunction and the mean in Figure 6. 

In Figure 7 we show the variance between the functions for 
a given basis set and order and the mean eigenfunctions 
for = 100. In general as the order increases the variance 
between the basis sets decreases. For some basis sets such 
as Chebyshev, Laguerre and Legendre polynomials there is 
already convergence even at a low order of ~ 40. At very 
low order A'^ < 10 the eigenfunctions are very noisy. Other 
basis sets such as Fourier and interpolation reach conver- 
gence at a high order ~ 100. The Fourier basis set is noisy 
since high order functions can introduce highly oscillatory 
modes into the eigenfunctions, this can also be seen in the 
A'^ = 60 panel in Figure 5. The tophat basis set requires an 

even higher order 1 1 . We note that there would be even better 
agreement between the non-tophat basis sets and the mean 
if we excluded the tophat basis set from this analysis. 



It can be seen that at Order = 40, for example, there is still 
a significant variation between basis sets. The variance be- 
tween basis sets decreases as the maximum order is increased 
but for the A^ = 4 function the tophat basis set varies sig- 



I We have not extended this calculation to even higher order 
since numerical effects involved in inverting (nearly singular) ma- 
trices of very high order start to become important, we choose a 
safe maximum order of 100 for this exercise. 
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4 CONCLUSION 

We have found there is some ambiguity in the use of eigen- 
decomposition in the literature to date and have aimed to 
clarify this situation. By investigating an arbitrary basis set 
expansion of w{z) we have found that the technique of find- 
ing the eigenfunctions of w{z) from the Fisher matrix is sen- 
sitive to the basis set choice and the order of the truncated 
expansion. 

We used the example of a weak lensing tomographic sur- 
vey to demonstrate how the eigenfunction decomposition is 
sensitive to these choices. We find that for a Euclid/DUNE- 
like survey the best constrained eigenfunctions could be de- 
termined to an accuracy of < 1% using lensing alone (no 
priors) . 

de Putter & Linder (2008) also presented a critique of the 
eigendecomposition approach which sheds further doubt on 
the technique's validity, especially at low truncated order. 
They note that whilst the noise term on an eigenfunction 
can be predicted it is in fact the signal-to-noise that will be 
important for future surveys. We note then that a low order 
expansion in a particular basis set is not robust and also has 
uncertain relevance in gauging the predicted performance of 
a future survey. 

Because of these ambiguities in the reconstruction of ■w{z) 
we recommend a move towards using physically motivated 
functional forms based on models. Physical parameters can 
be constrained and models compared using evidence calcu- 
lations (e.g. Trotta, 2007; Heavens et al., 2007 for a link with 
Fisher matrices). 

We summarise our conclusions here 



• The parameterisation of w{z) is an arbitrary basis set 
choice, and since we do not know the nature of dark energy 
in general any choice is valid. 

• When using orthogonal basis sets it is important to self- 
consistently include the weight function. 

• When one refers to "binning" one is actually referring 
to a tophat basis set. Interpolation is also valid but should 
not be confused with binning. 

• The set of eigenfunctions for a finite truncated basis set 
is not unique. There will always exists some eigenfunctions 
that are not common to all (truncated) eigenfunction sets. 

• When different basis sets are truncated to high order 
the best constrained eigenfunctions will tend to agreement. 

• When marginalising over extra cosmological parameters 
we show that finding the eigenfunctions of the full Fisher ma- 
trix leads to mixed non-'w{z) and w{z) eigenfunctions that 
are dependent on the non-i«(2;) parameterisation. We pro- 
pose an alternative in which the marginalisation affects the 
'w{z) eigenvalues only, but not the functions themselves, by 
finding the eigenfunctions of the w{z) sub-matrix. 

• To find agreement between the top ~ 4 eigenfunctions 
an order of at least 100 is required, for some basis set choices 
this is higher. 



The convergence of the eigenfunctions depends on the orig- 
inal choice of basis sets - some converge more quickly than 
others. We find that Chebyshev, Legendre and Laguerre 
polynomials reach convergence more quickly than interpola- 
tion, tophat or Fourier basis sets. 

Even with an order of 100, the residuals between the recon- 
structed eigenfunctions cause the associated errors to vary 
by as much as a factor of 1 - 10 for the best ~ 4 constrained 
functions. 

We note that the dark energy Figure of Merit suggested 
by Albrecht et al. (2009) uses a tophat basis set with 36 
parameters. 

We recommend that if an eigenfunction Fisher matrix anal- 
ysis is needed then a robust procedure should be followed 
with respect to marginalisation such that results are not af- 
fected by non-dark energy assumptions and priors and that 
a very high order of expansion is needed. The convergence 
can be tested by investigating a variety of basis sets - using 
only a single basis set is insufficient. 

The code used in the article will form part of the iCosmo 
open source package and will be available in version 1.2 at 
http : //www. icosmo . org. 
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APPENDIX A: EQUIVALENCE OF EIGEN 
AND JACOBIAN MATRICES 



In this Appendix we will show that the eigenmatrix is equal 
to a Jacobian matrix. The eigenvalue decomposition of a 
matrix can be written like 



QAQ' = F. 



(Al) 



The Jacobian transform of a Fisher matrix maps from one 
parameter set ai to another bi 



JAJ' =F 



(A2) 



where the matrix J in this case maps from the set 6 to a. 
Where A and F are defined as 



X 

= E 



da 



dhi 



(A3) 



and the signal is parameterised in as 
C{(j);x) = y^^ai(j}i{x) 

i 

C7(t/.;a;) = ^&iV.W (A4) 

i 

in F and A respectively. The elements of the Jacobian are 

(A5) 



When transforming from one basis set to another we can 
write the following 



1 



f{x)Mx)M{x)dx 



— ( VbmVm(a;) ) 0„(x)M(a;)dx 



(A6) 

i)mix)(t>nix)M^{x)dx (A7) 



so that the Jacobian can be written 
da„ _ 1 f 

where the weight M and constant c are for the (ji basis set 



If we assume that tp are eigenfunctions of the then they can 
be reconstructed using linear combinations of the original 
basis functions 



(A8) 



Since the decomposition of a matrix F to a diagonal matrix 
A in equations (Al) and equations (A2) is unique it can 
easily be seen that the eigenmatrix is in fact a Jacobian. 
This can be shown in the following way. The Jacobian is 
written like 



dan 
W 



1 



™ Cn J R 



■ip,n{x)(l)„{x)M^{x)dx. 



(A9) 



We construct new basis sets as linear combinations of the 
old so that 

dan _ 1 f 

this leads to 
dan 1 



ak(j>k{x)4)n{x)M4,{x)Ax 



(AlO) 



dbn 



Jmn 



k / (t>k{x)(l>n{x)M^{x)dx 



Cn J R 

T" Q mk^ knCn 
Cn 

Qmn- 



(All) 



APPENDIX B: NON-UNIQUENESS OF THE 
EIGENVALUES 



Here we will show that if the residual between the recon- 
struction of a basis set using another is non-zero then the 
eigenvalues obtained from each set via an eigenvalue decom- 
position of the Fisher matrices are not equal. 

The rotation matrix from basis (f) to tp can be written 

Jij = 4 / dxMx)^,{x)M^ix). (Bl) 

J R 

We now assume that the new basis functions can be writ- 
ten as the sum of some linear combination of the old basis 
functions plus some residual (equation 18) 



J., = - / dx(JtkM^) + Ri{x))cPj{x)M^ix) (B2) 



where the rotation matrix J rotates the old Fisher matrix 
to the new Fisher matrix assuming no residual. This can 
now be written like 



T 



4 / dxR,{x)(p,{x)M\x). 

J R 



(B3) 



In the following we will use subscripts for the basis for clar- 
ity. We now rotate from tp to %p then diagonalise the Fisher 
matrix of tp to get A^ this can be written as 
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We can also diagonalise F^^ directly to get 

If we assume that = A^ = A then comparing equations 
(B4) and (B5) the following must be true 

{J,f>-il> + T^^)J^A — J41A 

since we have assumed A^ = A^ the first pair of rota- 
tion matrices commute (however in general this is not true) 
J<j>iiJi,A ~ J 4, A and we have 

J4>A + T^4:J,pA = J,i>A (B7) 

which leads to 

7^1/. J1/.A = 

T^^ = 0. (B8) 

So if A^ = A^ then T^^ = 0. By deduction we can then 
assert that if T^^ 7^ then A^ 7^ A^ and the eigenvalues are 
not unique. 
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